Magnetism in Disordered Graphene and Irradiated Graphite 
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The magnetic properties of disordered graphene and irradiated graphite are systematically studied 
using a combination of mean-field Hubbard model and first principles calculations. By considering 
large-scale disordered models of graphene I conclude that only single-atom defects can induce fer- 
romagnetism in graphene-based materials. Preserved stacking order of graphene layers is shown to 
be another necessary condition for achieving a finite net magnetic moment of irradiated graphite. 
Ab initio calculations of hydrogen binding and diffusion and of interstitial-vacancy recombination 
further confirm the crucial role of stacking order in 7r-electron ferromagnetism of proton bombarded 
graphite. 

PACS numbers: 61.72.Ji, 61.80.Az, 75.75.+a, 81.05. Uw 



Graphene and related nanostructured materials are 
recognized as possible building blocks for prospective 
electronics technologies Q among which spintronics Q 
occupy a special place. Extraordinary magnetic proper- 
ties of carbon nanostructures resulting from their reduced 
dimensions have been predicted theoretically [1, 0> @j- 
This progress allowed to think of novel devices realiz- 
ing spintronics in practice [fjj. Although most of the 
predicted properties are still awaiting their experimen- 
tal confirmation, a discovery of room-temperature ferro- 
magnetism in proton irradiated graphite [7[ continues to 
stimulate research in the field of magnetic carbon. 

Recent experimental investigations have revealed that 
the magnetic order in proton bombarded graphite has 
two-dimensional (i.e. graphene-like) character 0] and 
originates from the carbon 7r-electron system rather than 
from intrinsic or introduced impurities These re- 



sults are supported by theoretical calculations [lfj, [ll| 
which show the crucial role of defect-induced quasilocal- 
ized states 12, Il3| in producing a magnetic order sta- 



ble at high temperatures. Although magnetic properties 
of defects in graphene have been widely studied by first 
principles methods, the use of different periodic mod- 
els of limited size often resulted in varying outcomes 

Q EE 0, G3- This calls for a stud y employing 
realistic disordered models. The origin of ferromagnetic 
order in graphene which generally favors antiferromag- 
netic coupling between the neighboring atoms 0, as 
well as the role of various possible defects in inducing 
magnetism remain elusive. 

In this Letter, I address magnetic properties of disor- 
dered graphene and irradiated graphite by using a combi- 
nation of mean-field Hubbard model and first principles 
calculations. The Hubbard model calculations of realis- 
tic models of disordered graphene show that only single- 
atom defects such as vacancies or chemical functionaliza- 
tions unequally distributed over the two sublattices may 
result in a net magnetic moment. Further first principles 



calculations point out that the stacking order of graphite 
lies at the origin of ferromagnetism in proton bombarded 
graphite. 

To study the magnetic properties of disordered 
graphene I employ the mean-field approximation of the 
Hubbard model which proved to describe carbon tt- 
electron systems in good agreement with first principles 
calculations [B| . The corresponding Hamiltonian reads 



H 
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h.c] + U^2n it7 (ni^ a ), (1) 



where the first part is the single-orbital tight-binding 
Hamiltonian while the second part accounts for the on- 
site Coulomb repulsion. In this expression ci a (c| CT ) an- 
nihilates (creates) an electron with spin a at site i and 
stands for the pairs of bonded atoms. The expec- 
tation values of the spin-resolved density rii a = c\ a Ci a 
on atom i are obtained from the eigenvectors of Ji. The 
self-consistent solution of the problem provides the spin 
densities Mi = (n^ — «ij,)/2 on each atom i. 

The models of disordered graphene arc generated by 
distributing randomly point defects in a large 20 x 20 
supercell (800 atoms). The calculations are performed 
for 64 different random distributions in order to pro- 
vide reliable statistical sampling. This computational 
protocol yields the mean magnetic moments converged 
within ~1% both with respect to the supercell size and 
the number of configurations. Similar approach was suc- 
cessfully applied to study vacancy-induced magnetism 
in oxide materials (l8| . I consider the range of de- 
fect concentration x = N^/N = 0.01 — 0.1 with 
out of N atoms or bonds being affected. The values 
U/t = 1.0,1.33,1.67,2.0 are studied in order to inves- 
tigate the role of on-site Coulomb repulsion magnitude. 
While the value of hopping integral t « 2.7 eV is well 
established for graphene, there is a growing debate re- 
garding the on-site Coulomb repulsion U. It is worth 
pointing out that the magnitude of U was inferred from 
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magnetic resonance studies of neutral soliton states in 
irans-polyacetylene a one-dimensional bipartite sp 2 
carbon system closely related to graphene. The esti- 
mated U = 3.0- 3.5 eV {U/t = 1.1 - 1.3) [2fj] closely 
corresponds to U/t = 1.33 considered here. 

The empirical value of U/t can also be compared to the 
values obtained by establishing a relation between the 
results of first principles and mean-field Hubbard model 
calculations [21( . The local density approximation (LDA) 
leads to U/t » 0.9 while the generalized gradient approx- 
imation (GGA) gives U/t « 1.3 in good agreement with 
the empirical value. It is somewhat cautionary that the 
results provided by generally more accurate hybrid den- 
sity functionals correspond to U/t « 2.0 which is very 
close to the critical value U/t « 2.2 at which graphene ac- 
cepts an antiferromagnetic ground state [22} ■ Discrepan- 
cies between the results of hybrid density functional cal- 
culations and experimental data for irans-polyacetylene 
have also been reported j23j. This urges to reconsider 
performance of these methods for describing magnetic 
sp 2 carbon systems. 

Disorder in graphene induced by chemical treatment 
or irradiation with ions or electrons can be conveniently 
classified by the following three types of point defects. 
A carbon atom removed from the sp 2 lattice (iW'va- 

or 
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cancy') corresponds to either true vacancy 1 
the rehybridization into the sp 3 configuration. Such re- 
hybridization may result from chemical functionalization 
(e.g. hydrogen chemisorption 1(| 3]) or binding of the 
interstitial carbon atoms in irradiated graphite 2lL 26j |. 
A multiple-atom 'vacancy' can be viewed as an ensem- 
ble of neighboring single-atom 'vacancies' . In some cases 
a bivalent chemical functionalization results in breaking 
the bond between the pairs of adjacent atoms without 
their rehybridization [27]]. This leads to so-called bond 
dilution (ii). While (i) and (ii) give rise to undercoordi- 
nated carbon atoms, topological defects (iii) such as the 
Stone- Wales defect [28( maintain the coordination num- 
ber. This type of disorder can be described as a permu- 
tation of connectivity involving two or more neighboring 
atoms. However, no magnetic moments due to the pres- 
ence of the Stone- Wales defects were observed in this 
study as well as in previous ab initio calculations [Hj]. 
Below, I focus only on the magnetism induced by the 
'vacancies' and the bond dilution disorder. 

The Hubbard model solutions for all configurations 
considered are characterized by the antiparallel orien- 
tation of magnetic moments on the two sublattices. 
Figs. QJa,b) show the sublattice-resolved average mag- 
netic moments (M?) (77 = A, B) induced by varying 
concentrations of 'vacancies' for different values of U/t. 
For the case of defects equally distributed over the two 
sublattices [Fig.QJa)] (M^) and (Mf) compensate each 
other resulting in an antiferromagnetic state with zero 
total magnetic moment M = N((M^) + (Mf))/2 = 0. 
This result can be understood within the Lieb theorem 




FIG. 1: (Color online) Average magnetic moments of carbon 
atoms in A and B sublattices of graphene as a function of 'va- 
cancy' concentration for different U/t values. The 'vacancies' 
are either (a) distributed equally among sublattices or (b) be- 
long to sublattice B only, (c) Spin density distribution in a 
representative configuration with 'vacancy' defects (x = 0.03) 
in sublattices A (A) and B (T). 



[291 ] which relates the total magnetic moment of a half- 
filled bipartite system to the difference of number of 
atoms in the two sublattices. In our case the two sub- 
lattices contain equal number of atoms resulting in zero 
net magnetic moment. The sublattice-resolved magnetic 
moments (M?) a i 1 (0 < 7 < 1) show sub- linear de- 
pendence on x and tend to increase with increasing the 
magnitude of U/t. In particular, for U/t = 1.33 the 
dependence on 'vacancy' concentration x can be accu- 
rately fitted by \(M?)\ w 0.234 • x - 78 . For the case of 
'vacancies' in sublattice B only the total magnetic mo- 
ment per unit cell M — x [open circles in Fig. QJb)]. 
This is again in strict agreement with the Lieb theo- 
rem predicting 2M = — Nq. The magnetic state 
originates from the Hund's rule population of — N B 
zero-energy defect-induced quasilocalized states in sub- 
lattice A. However, these singly-occupied states induce 
exchange polarization of the fully-occupied ones leading 
to a significant antiparallel magnetization of the atoms 
in sublattice B. The induced magnetic moments follow 
an approximate dependence (Mf) w —0.159 • x 069 . 

Figure[TJc) illustrates the distribution of magnetic mo- 
ments induced by the 'vacancy' disorder in a representa- 
tive configuration (x = 0.03) with defects in both sublat- 
tices. Magnetic moments develop in the nanometer-size 
regions characterized by the local domination of 'vacan- 
cies' in one of the sublattices [1 and 2 in Fig. HJc)]. In 
these regions one can trace distinctive triangular patterns 
associated with individual quasilocalized states [l2| . The 




FIG. 2: (Color online) (a) Average magnetic moments of car- 
bon atoms in A and B sublattices of graphene as a function 
of bond dilution for different U/t values and (b) spin den- 
sity distribution in a representative configuration (x = 0.10; 
U/t= 1.33). 



orientation of magnetic moments is determined by the 
antiferromagnetic coupling between the two sublattices. 
In contrast, the proximity of two defects pertaining to op- 
posite sublattices [3 in Fig. [ljc)] produces no local mag- 
netic moments. This can be explained by the fact that 
the quasilocalized states on different sublattices form a 
pair of bonding/antibonding states [Til ] . and the energy 
splitting increase with decreasing the distance between 
the two 'vacancies'. Such dispersion provides a mecha- 
nism for escaping the low-energy instability which is pre- 
ferred to the spin polarization. Increased contribution of 
this channel with increasing defect concentration explains 
the sublinear dependence of induced magnetic moments 
on x for the case of 'vacancies' in both sublattices. Both 
mechanisms contribute to stabilizing low-energy states. 
Indeed, band gaps increasing from 0.14 eV to 0.29 eV 
are found in the range of concentrations x = 0.01 — 0.1 
and U/t = 1.33. This agrees with the recent experimen- 
tal observation of insulating behavior in proton irradiated 
graphite [HJ. 

In contrast to the 'vacancy' disorder, bond dilution 
does not change the number of atoms and, thus, can not 
produce a net magnetic moment. Nevertheless, antifer- 
romagnetic correlations develop as a result of bond dilu- 
tion [Fig.|2fa)]. At low values of bond dilution (x < 0.05 
for U/t = 1.33) only very small average magnetic mo- 
ments develop as a results of occasional occurrence of 
highly undercoordinated atoms. For larger values of x 
spin polarization occurs in extended regions of higher lo- 
cal concentration of broken bonds as shown in Fig. [2](b) 
for x = 0.1 and U/t = 1.33. Such cumulative behavior 
is similar to the quantum phase transition predicted for 
hexagonal graphene nanoislands [f|. 

One can now conclude that a macroscopic magnetic 
moment may appear only in the presence of 'vacancy' 
disorder unequally distributed over the two graphene 
sublattices. In bulk graphite and multilayer graphene 
the sublattices are inequivalent due to the stacking or- 
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FIG. 3: (a) Inequivalent carbon atoms (A and B) in graphite, 
(b) Possible pathways for the diffusion of chemisorbed hydro- 
gen in the ab plane of graphite, (c) Relative energies of the 
potential energy surface minima and transition states for the 
diffusion of chemisorbed hydrogen in graphite. 



der [Fig. [Sfa)]. Irradiation of graphite with high-energy 
ions initially produces (i) vacancy-interstitial pairs due 
to knock-on recoils and (ii) chemisorbed ions stopped by 
the bulk material. It was found that proton bombard- 
ment of graphite leads to much stronger magnetic sig- 
nal than the He-ion irradiation where no chemisorption 
takes place. This suggests that the hydrogen chemisor- 
tion plays an important role in the onset of magnetism 
in proton bombarded graphite. 

The effect of stacking order in graphite is further inves- 
tigated by means of first principles calculations. I consid- 
ered a 4 x 4 x 1 supercell of ABA stacking graphite within 

" M. 



the GGA density functional theory framework [31 



Chemisorption of a proton in position B was found to be 
0.16 eV lower in energy than in position A [Fig. [3£b,c)]. 
This is explained by the steric repulsion between the pro- 
ton chemisorbed on atom A and carbon atom A' of the 
adjacent graphene layer. A magnetic moment of 1 hb lo- 
calized in the functionalized graphene sheet is observed 
for the hydrogen chemisorption in both inequivalent po- 
sitions of the crystalline lattice of graphite. This is in full 
agreement with the previous ab initio calculations T3, 12 1 
and the present Hubbard model results for graphene. 
The barrier for the migration of proton initially trapped 
in position A to the neighboring position B is 0.88 eV 
(path A — * B in Fig. GO^c)). Such migration can take 
place at the experimental conditions. The overall barrier 
for the in-plane diffusion of protons in graphite (paths 
A -> B, B -> B' and A' -> B in Fig. E£b,c)) defined by 
the lowest barrier of path A — > B is 1.04 eV. Diffusion 
in the direction of c-axis requires much higher activation 
energy of ~4 eV. Free diffusion of protons in graphite will 
result in either saturation of vacancy dangling bonds 14 1 
or in the formation of energetically favorable aggregates 
with both sublattices being populated 17[. Magnetic 



moments of the hydrogen aggregates quench due to the 
large splitting of quasilocalized states in the complemen- 
tary sublattices. In addition, the condition of preserving 
the stacking order of graphite limits the maximum irradi- 
ation dose and intensity. Subsequently, there exist only 
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relatively narrow range of conditions at the sublattice- 
discriminating hydrogen chemisorption can be achieved. 
This accords with the fact magnetic moments are induced 
only in the narrow ring surrounding the irradiated spot 

Similarly to hydrogen chemisorption, single-atom va- 
cancies produced by the knock-on recoils induce low- 
energy quasilocalized states [loj |. In contrast, interstitial 
carbon atoms in the energetically preferred configuration 



331 ] lead to the rehybridization of the pairs of neigh- 
boring atoms in the adjacent graphene layers. This situa- 
tion is equivalent to the case of di- vacancies providing no 
low-energy states required for developing magnetic mo- 
ments. Cross-sections for the momentum transfer due 
to knock-on collisions are likely to be equal for both A 
and B carbon atoms in graphite. However, the stack- 
ing order of graphite has strong influence on the kinetics 
of interstitial-vacancy recombination. Instantaneous re- 
combination of low-energy recoils was found to be signif- 
icantly more likely for the atoms in position A (2||. At 
longer time scales, diffusing interstitials encounter differ- 
ent barriers for the recombination with vacancies in these 
two positions. First principles calculations show different 
minimum transition barriers of 0.04 eV and 0.22 eV for A 
and B vacancies, respectively. Both values are below the 
diffusion barrier of interstitial atoms (~0.5 eV [34]) sug- 
gesting that only very small difference in populations of 
the two sublattices by vacancy defects can be achieved in 
practice. Nevertheless, both vacancies and chemisorbed 
hydrogen atoms are found to populate preferably sublat- 
tice B thus producing fcrrimagnetic order with the spin 
polarization larger for sublattice A. 

In conclusion, I have shown the crucial role of 
single-atom defects in combination with a sublattice- 
discriminating mechanism for developing ferromagnetic 
order in graphene-based materials. In graphite the role 
of such mechanism is played by the stacking order of 
graphene layers. This suggestion is confirmed by the ab 
initio investigation of hydrogen chemisorption and va- 
cancy defects in proton bombarded graphite. The re- 
sults of this study pave a way for tailoring carbon based 
magnetic materials by means of irradiation and chemical 
treatment of graphite and other graphene derivatives. 
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